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model  calculations. 
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SUMMARY 


Orthophosphoric  acid  Is  commonly  resorted  to  when  physical  property  data 
are  required  to  model  phosphorus  smoke.  It  has  been  shown  previously  that  optical 
constants  of  orthophosphoric  acid  solutions  provide  reasonably  accurate 
approximations  to  the  extinction  properties  of  phosphorus  smoke  in  the  visible,  near 
infrared  ( IR) ,  and  3-  to  5-p®  IR  region,  but  not  in  the  8-  to  12-um  region. 
Furthermore,  studies  on  the  high-temperature  reaction  between  phosphorus  pentoxide 
and  water  vapor,  a  precursor  reaction  to  the  formation  of  a  phosphorus  smoke, 
indicates  the  chemical  composition  of  the  smoke  droplets  is  a  mixture  of  condensed 
phosphoric  acids.  From  these  results,  it  is  clear  that  phosphorus  smoke  is  not 
composed  entirely  of  orthophosphoric  acid. 

It  is  the  intent  of  this  study  to  delineate  the  expected  range  of  errors 
resulting  from  the  employment  of  an  orthophosphoric  acid  thermodynamic  model  for  the 
prediction  of  the  equilibrium  properties  of  a  phosphorus  smoke  under  tactical 
conditions.  Subsequently,  the  thermodynamic  properties  of  the  solution  droplets 
comprising  a  phosphorus  smoke  are  formulated  by  accounting  for  the  condensed  nature 
of  the  phosphoric  acids.  By  considering  the  condensed  phosphoric  acids  as  a 
homologous  series,  the  solution  thermodynamics  of  the  condensed  phosphoric  acids  is 
aproximated  from  the  orthophosphoric  acid  water  activity  data  by  applying  a  solute 
number  reduction  scheme. 

Asserting  the  acid  mixture  is  a  combination  of  pyrophosphor ic  acid  and 
orthophosphoric  acid,  a  mixture  model  is  proposed  for  the  representation  of  a 
phosphorus  smoke.  Comparison  of  mixture  model  data  and  experimental  data  indicates 
the  mixture  model  predicts  the  acid  concentrations  and  yield  factors  of  a  phosphorus 
smoke  within  5%.  This  represents  a  significant  improvement  over  the  25%  errors 
resulting  from  the  orthophosphoric  acid  model  predictions. 


PREFACE 


The  work  described  In  this  report  was  authorized  under  Project  IT  161 102 A7  1  A, 
C/B  Defense  and  General  Investigations;  Scientific  Area  5,  Aerosol/ Obscuration  Science.  This 
work  was  conducted  fiotn  January  to  November  1979.  The  background  data  are  contained  in 
notebook  9990. 

Reproduction  of  this  document  In  whole  or  in  part  is  prohibited  except  with 
permission  of  the  Com  m ander/ Director,  Chemical  System o  Laboratory,  ATTN:  DRDAR-CLJ-R, 
Aberdeen  Proving  Ground,  Maryland  21010;  however,  the  Defense  Technical  Inform  ation  Center 
and  the  National  Technical  Information  Service  are  authorized  to  reproduce  the  document  for 
United  States  Government  purposes. 


CONTENTS 


Pago 


]  INTRODUCTION  .  9 

2  THERMODYNAMIC  PROPERTIES  OF  ORTHOPHOSPHORIC  ACID  DROPLETS .  9 

3  REFORMULATION  OF  PHOSPHORUS  SMOKE  THERMODYNAMIC  PROPERTIES .  14 

3.1  Solution  Thermodynamics  of  Condensed  phosphoric  Acids  .  13 

3.2  Equivalent  H^PO^  phosphorus  Smoke  Acid  Concentration  . .  17 

3.3  Real  Phosphorus  Smoke  Acid  Concentration  . .  19 

3.4  Yield  Factor  of  a  Phosphorus  Smoke  .  21 

4  CONCLUSIONS  .  22 

LITERATURE  CITED  .  23 

GLOSSARY  .  25 

DISTRIBUTION  LIST  .  27 


LIST  OF  FIGURES 

Figure 

1  Acid  Concentrations  of  a  Phosphorus  Smoke  Versus  Relative 

Humid  i  ty  . . .  12 

2  Yield  Factors  of  a  Phosphorus  Smoke  Versus  Relative  Humidity.. .  14 

3  Water  Activity  of  Various  Condensed  Phosphoric  Acid  Solutions 

in  a  Function  of  H-jPO^  Association  .  16 

4  Equivalent  H^PO^  Acid  Concentation  of  a  Phosphorus  Smoke 

Versus  Degree  of  Acid  Association  . . .  18 

5  Percent  Deviation  hetween  Equivalent  H  PO^  and  Real  Acid  Concen¬ 

tration  as  a  Function  of  AcL  Association.. .  20 

TABLE 

Dependence  of  Kelvin  Effect  on  Droplet  Diameter  for  Sized 

Representative  of  a  Phosphorus  Smoke .  10 


7A 


AN  IMPROVED  THERMODYNAMIC  MODEL  FOR  PHOSPHORUS  SMOKE 


1.  INTRODUCTION 

To  accurately  predict  the  condensational  growth  of  a  phosphorus  smoke  and 
its  attendant  obscuration  efficiency,  detailed  knowledge  of  the  solution 
thermodynamics  of  the  acid  droplets  is  required.  Previous  attempts  to  model  the 
thermodynamic  properties  of  a  phosphorus  smoke ^  have  been  based  on  the  assumption 

that  the  generic  acid  of  the  solution  droplets  is  orthophosphori.  acid  (H  PO  ). 

9  3  4 

However,  recent  transmission  measurements  by  Milham  etal.  have  Indicated  the 
optical  properties  of  a  phosphorus  smoke  are  distinctly  different  from  these  of  a 
comparable  H^PO^  mist.  These  findings  coupled  with  the  high-teraperature  phosphorus 
pentoxide  (P4010)  hydration  studies  by  Kapishe  3  .  learlv  indicate  that  a  phosphorus 

smoke  is  not  composed  entirelv  of  H,P0, . 

'  3  4 

It  is  the  purpose  of  this  study  to  address  the  expected  range  of  errors 
resulting  from  the  representation  of  a  phosphorus  smoke  by  a  thermodynamic  model 
that  is  based  on  the  solution  propert.es  of  H^PO^ .  Having  demonstrated  the 
inadequacies  of  the  previous  model,  a  more  realistic  thermodynmami c  model  is 
formulated  by  accounting  for  the  condensed  nature  of  the  droplet  phosphoric  acids. 
Based  on  the  high  temperature  P^Oj  hydration  studies  by  Kapishev,  it  is  postulated 
the  acid  components  of  a  phosphorus  smoke  are  75%  py rophosphoric  acid  and  25%  !L  PC, 
by  weight.  The  solution  thermodynamics  of  an  arbitrary  mixture  of  condensed 
phosphoric  acld6  is  established  from  the  H^PO^  water  activity  data  through  a  solute 
number  reduction  scheme.  Comparison  of  mixture  model  synthetic  data  and  H,PO,  mode! 
predictions  with  experimental  data  taken  from  a  phosphorus  smoke  shows  that., 
although  errors  of  25%  are  associated  with  tne  H^PO^  model,  maximum  errors  of  only 
5%  result  from  the  mixture  model  predictions.  The  theoretical  development  of  a 
thermodynamic  model  for  the  condensed  phosphoric  acids  is  believed  to  be  a  more 
robust  representation  of  phosphorus  smokes  than  previous  empirical  attempts. 

2.  THERMODYNAMIC  PROPERTIES  OF  ORTHOPHOSPHORIC  ACID  DROPLETS 

The  acid  droplets  are  proposed  to  be  mixtures  of  volatile  water  and 
nonvolatile  acids.  In  general,  the  water  partial  pressure  over  the  droplet  mi; face 
is  a  function  of  the  droplet  diameter  and  droplet  acid  concentration.  At  os  ant 
temperature  T  and  total  pressure  P,  the  droplet  saturation  ratio  S  is  found  by 
minimizing  the  solution  droplet  Gibbs  free  energy.  The  droplet  saturation  ratio, 
the  degree  of  water  saturation  at  the  droplet  surface,  is  given  by 
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The  first  term  on  the  right-hand  side  is  referred  to  as  the  Kelvin  effect,  a  vapor 
pressure  elevation  effect.  The  Kelvin  effect  causes  verv  small  droplets  to  have 
extremely  high  vapor  pressures,  whereas  the  larger  droplets  are  affected  to  only  a 
small  extent.  The  table  illustrates  the  increase  in  vapor  pressure  over  phosphoric 
acid  droplets  due  to  its  increasing  curvature.  It  is  evident  that,  for  droplets 
with  diameter  greater  than  0.3  pm,  the  Kelvin  effect  represents  less  than  a  1° 
elevation  in  the  droplet  vapor  pressure. 


Table.  Dependence  of  Kelvin  Effect  on  Droplet  Diameter 
for  Sizes  Representative  of  a  Phosphorus  Smoke 


Exp  Aov^/RTd 

d 

um 

1.23 

0.01 

1  .02 

0.10 

1.01 

0.30 

1.00 

1.00 

Cascade  inpactor  measurements  by  Gillespie  and  Johnstone'  show  that  only  those 
droplets  with  diameter  greater  than  0.3  um  contribute  significantly  (>99%)  to  the 
phosphorus  smoke  mass  distribution.  Consequently,  thermodynamic  properties  of  the 
non-Kelvin  phosphoric  acid  droplets  are  completely  specified  by  the  relation 

S*at  (2) 


The  water  activity  a^  of  the  acid  solutions  is  defined  as  the  ratio  of  the  water 
pressure  over  the  acid  solution  to  the  saturation  water  vapor  pressure  at  the 
solution  temperature.  The  nonvolatile  acid  depresses  the  water  vapor  pressure  and 
the  droplet  saturation  ratio  will  be  less  than  unity. 

From  water  vapor  pressure  data  over  a  discrete  set  of  H,P0^  solutions  as 
given  by  Mellor,5  the  water  activity  of  H^PO^  solutions  can  be  analytically  related 
to  the  water  mole  fraction.  Toward  the  aim  of  finding  a  continuous  analytic 
relationship,  a  nonlinear  regression  analysis  is  applied  to  the  water  activity 
data.  Within  a  2%  error,  the  water  activity  of  H^PO^  solutions  can  be  represented 
by  the  multivalued  function 


3  2 

n]  =  -7 .81 2Xj  +  1 8 . 828Xj  -  1 2 . 70  3x  t  +  2.711 

a  j  =  -6.343Xj  +  13.649xl  -  6.308 


0.10  <  a 


<0.75 


1  0 


0.75  <a  j  <0.99 


(3) 


It  has  been  show  previously  by  Rubel^  that  the  acid  nuclei  produced  by 
the  combustion  of  elemental  phosphorus  condensatlonally  grow  to  a  Btable  equilibrium 
within  seconds  of  formation.  These  calculations  have  been  experimentally  supported 
by  the  acid  concentration-residence  time  measurements  of  Gillespie  and  Johnstone. ^ 
Stable  equilibrium  for  a  solution  droplet  is  defined  as  the  condition  where  the 
droplet  saturation  ratio  equals  the  environmental  relative  humidity.  Therefore,  the 
acid  droplets  have  achieved  diffuslonal  equilibrium  under  tactical  conditions  and 


By  inverting  equation  4,  the  water  mole  fraction  can  be  defined  in  terms  of  t he 
relative  humidity  ¥  as 


X  »  0.632  cos  [cos  (1.175-1.928*0/3  +  240]  +  0.803  O.IK'KO.75 


2.152  -  (0.655  -  0.632Q 


1/2 


(4) 


0.75<¥ <0.Q9 


fraction  x 


From  the  defining  relationship  between  acid  fraction  F  and  water  mole 
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= _ 1.000  -  F 

(1  *  1.00  -  0.816F 


(5) 


the 


acid  fraction  of  an  H^PO^  droplet  can  be  given  as 


0.197  -  0.632C 
0.345  -  0.516C 
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and 


F  - 


-0.076  +  0.500  (0.655  -  0.632*0 
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It  is  noteworthy  that  the  duplet  acid  concentration  F  is  independent  of 
the  droplet  diameter  and  is  solely  a  Function  of  relative  humidity.  As  a 

consequence,  all  acid  droplets  are  characterized  by  the  same  acid  concentration  at  a 
given  relative  humidity. 

Anderson*  measured  the  acid  «.  oncentrat ion  of  phosphorus  smoke  liquid  from 
20°;  to  90%  relative  humidity.  The  measurement  involved  titrat'ng  the  liquid  sample 
against  a  standard  0.01  N  NaOH  solution.  The  H^PO^  molecule  strongly  dissociates 
one  hvdrogen  ion  allowing  tor  a  substitution  bv  the  sodium  ionr  Since  H,P0,  i 

3  a 

primarily  monoprotic ,  each  mole  of  Ns  OH  dissolved  is  equated  to  a  mole  of  H^PC,.  • 

From  the  total  mass  of  sample  and  the  moles  of  H-PO,  present,  the  H.PO 

.3  ^  3 

concentration  is  determined. 

Figure  1  shows  the  differences  between  measured  acid  fractions  and  those 
values  predietd  by  the  H^PO  model.  The  ortho  model  consistently  underestimates  the 
acid  concentrations  over  the  relative  humidity  range  from  20%  to  90%.  F.rrors  vary 
from  14%  at  25%  relative  humidi  _y  to  a:  high  as  25%  at  90%  relative  humidity.  This 
would  suggest  t ha  ortho  model  overestimates  the  water  absorbing  characteristics  of 
the  acid  droplets.  Since  the  specific  solution  acid  controls  the  degree  of 
condensational  growth,  the  disparities  indicate  the  generic  acid  is  not  H^PO^. 


py 


Fice.re  1.  Acid  Concentrations  of  a  Phosphorus  Smoke 
Versus  Relative  Humidity 


*  Anderson,  p.  H .  Chemical  Systems  Laboratory,  Research  Division,  Aberdeen  proving 
St  -end,  Maryland..  Private  communication.  1978. 
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A  further  indication  of  the  inadequacies  of  v'e  ortho  model  Is  the 
significant  differences  between  measured  and  predicted  yield  factors.  The  yield 
factor  (y.f.)  is  defined  as  the  ratio  of  the  final  aerosol  mass  to  the  mass  of 
starting  smoke  generating  material.  If  the  omoke  material  is  pure  phosphorus  and 
each  phosphorus  atom  in  the  starting  material  is  aerosolized,  then  the  phosphorus 
smoke  yield  factor  is^ 

M(H  PO  ) 

y,f *  M(P)*  F  (7} 

where  MfH-PC^)  is  the  molecular  weight  of  the  H^PO^  molecule  and  M(P)  is  the  mass  of 
the  phosphorus  atom. 

Fquation  9  illustrates  the  mass  increase  of  a  starting  phosphorus  atom  as 
it  is  chemically  converted  to  H^PO^  which  then  experiences  water  absorption 
according  to  its  thermodynamic  properties.  To  predict  the  mass  of  phosphorus  smoke- 
generated  by  a  pure  phosphorus  starting  material,  the  initial  mass  is  scaled  upward 
by  the  appropriate  relative  humidity  dependent  yield  factor.  In  general,  all  of  the 
starting  phosphorus  is  not  aerosolized  due  to  perhaps  burning  efficiencies 
associated  with  the  design  of  the  starting  material.  The  phosphorus  aerosol  i  zn  t  i  on 
parameter  C,  which  is  normally  less  than  unity,  is  defined  as  the  ratio  of 
phosphorus  in  the  aerosol  state  to  the  phosphorus  in  the  s  arting  material. 
Furthermore,  a  particular  smoke  material  may  be  composed  of  phosphorus  and 
nonphosp'uorus  material.  Consequently,  the  realized  yield  factor  is  expressed  in 
terms  of  the  material  yield  factor  as 

y.f.  (realized)  =  n  *  C  *  y.f.  (material)  (8) 


where  n  is  the  percent  mass  of  phosphorus  in  the  starting  material.  Tarnovo' 
measured  the  material  yield  factor  of  a  phosporus  smoke  for  relative  humidities 
beLwet-n  *.()%  and  70%  relative  humidity.  Figure  2  depicts  the  disparities  between 
measured  and  ortho  model  predicted  yield  factors.  The  ortho  model  overestimates  the 
measured  yield  factors  by  as  much  as  25%.  The  overestimation  of  mass  is  a 
combination  of  generic  acid  and  attendant  solution  thermodynamics  misrepresentation. 

Since  the  disparities  are  significant,  it  may  be  concluded  that  the 
model  cannot  satisfactorily  describe  the  thermodynamic  properties  of  a  phosphorus 
smoke.  In  fact,  the  failure  of  the  ortho  model  is  not  totally  unexpected  if  oiu 
considers  the  high  temperature  hydration  of  P^Ojq,  a  precursory  reaction  to  t  in¬ 
formation  of  a  phosphorus  smoke.  Although  it  is  commonly  assumed  that  the  P.0 
hydration  product  is  H^P1-^ ,  Kapishev^  has  found  that,  for  temperatures  about  2!M'"o, 
the  hydration  reaction  results  In  a  mixture  of  condensed  phosphoric  acids.  onlv 
after  days  will  these  condensed  acids  revert  to 
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Therefore,  to  accurately  model  the  condensational  growth  of  a  phosphorus 
smoke,  thermodynamic  data  arc  required  on  the  solution  thermodynamics  of  the 
condensed  phosphoric  acids.  Unfortunately,  the  thermodynamic  data  on  the  condensed 
phosphoric  acids  are  sparse.  However,  the  water  activity  of  the  condensed 
phosphoric  acids  may  he  heur i at  leal ly  formulated  by  invoking  key  assumptions. 

1.  RFF0RMIM.ATION  OK  PHOSPHORUS  SMOKE  THERMODYNAMIC  PROPERTIES 

The  combustion  of  elemental  phosphorus  establishes  a  severe  super- 
saturation  of  P^O  q  vapor.  The  nucleation  of  the  oxide  is  most  probably  governed  by 
the  sublimation  characteristics  of  the  having  a  sublimation  temperature  of 
around  250°C.^  Once  nucleated,  the  oxide  chemisorbs  water  vapor  and  converts  to  a 
phosphoric  acid.  kaptshev^  experimentally  analyzed  the  high-temperature  hydration 
reaction  between  F'  Oj^  ar*c'  water  vaoor.  Employing  chromatographic  methods,  the 
researchers  concluded  that  the  resultant  liquid  was  a  mixture  of  condensed 
phosphoric  acids  whose  percent  composition  was  temperature  dependent.  For 
chemisorption  at  20  c  the  percent  weight  composition  of  the  mixture  was  determined 
to  be  25"!  H,?0^  and  7 5%  py rophosphor lc  acid.  This  acid  composition  will  represent 
the  hygroscopic  solute  employed  In  the  following  analysis  of  phosphorus  smoke 
growth. 
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Yield  Factors  ol  a  Phosphorus  Smoke  Versus  Relative  Humidity 


Solution  Thermodynamics  of  Condensed  Phosphoric  Acids 


t.I 


The  homologous  series  of  condensed  phosphoric  acids  can  be  represented  by 
Che  chemical  formula^ 


i(H3P04)  -  (i-i)H20 


(9) 


where  i  represents  the  order  of  phosphoric  acid  association,  i.e.,  i  «  1  ♦  H^PO^  and 
i  ■  2  +  orthophosphorlc  and  pyrophosphoric  acid,  respectively.  Therefore, 
the  condensed  phosphoric  acids  can  be  thought  of  as  physi ochemical  molecular 
associations  of  H^PO^  molecules,  except  for  small  differences  in  molecular  weight. 

The  water  activity  of  a  solution  is  a  strong  decreasing  function  of  the 
relative  number  of  solute  (acid)  molecules  in  solution.  Since  the  associated 
solution  Is  characterized  by  a  smaller  number  of  solute  molecules,  the  associated 
solution  water  activity  is  greater  than  the  unassociated  solution  water  activity  for 

Q 

the  same  mass  of  acid.  This  conclusion  was  numerically  supported  by  the 
calculations  of  Ogsten  and  Winzer,^  who  considered  the  effects  of  associating 
solutes  on  the  osmotic  pressure  of  a  solution.  Consequently,  for  the  same  mass  of 
acid,  the  associated  solution  will  absorb  a  smaller  amount  of  water  than  the 
unassociated  solution.  The  associated  acid  concentration  will  be  correspondingly 
greater  than  the  unassociated  acid  concentration.  Such  a  circumstance  would  explain 
the  disparities  depicted  figure  1. 


The  acid  mixture  solvent  mole  fraction  may  be  written  as 

/ 

n 

Xi 


H2° 


4  n'«2o  +  ini 


(10) 


where  j!n^  is  the  number  of  condensed  phosphoric  acids  resulting  from  the  association 
of  |inj  H3PO4  molecules.  This  follows  from  the  hydration  process  of  condensed 
phosphoric  acids  where  the  ith  order  condensed  phosphoric  acid  produces  1  molecules 
of  H3PO/,.  The  mixture  solvent  mole  fraction  can  be  expressed  in  terms  of  the  ortho 
model  parameters  by  the  relations 


h2c 


1^0 


and 


l1"! 


(ID 


where  the  association  parameter  a  represents  the  degree  of  H3PO4  molecule 
association.  Substitution  of  equation  11  into  equation  10  gives 


xi  ‘  * /u  +  PV"^ 


(12) 
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In  general,  deviations  from  Raoult's  law  result  from  solute-solvent  interactions.® 
Due  to  a  strong  hydrogen  association,  H^PO^  appears  electronegative  In  solution. 
Consequently,  high  dipole  moment  water  molecules  will  interact  with  ionized  V^PO^ 
molecules.  It  is  asserted  that  this  mode  of  solvent-solute  interaction  is  the  basis 
for  the  nonideality  of  H^PO^  solutions  and  governs  the  nonlinear  dependence  of  water 
activity  on  water  mole  fraction  (equation  3).  The  condensed  phosphoric  acids 
dissociate  similarly  and  are  assumed  to  interact  with  the  water  molecules 
analogously.  Therefore,  the  functional  dependence  of  the  acid  mixture  solvent  mole 
fraction  xj  is  specified  by  equation  3  in  this  study. 

Figure  3  graphically  illustrates  the  dependency  of  solution  water 
activity  on  mole  ratio  of  H^PO^  (before  association)  to  water  for  various  degrees 
of  acid  association.  The  degree  of  association  increases  with  the  association 
parameter  _a_,  so  that  *  1  is  a  noncondensed  pure  H^PO^  solution  and  •  4  is  a 
highly  condensed  tet ra-metaphosphor ic  acid  solution.  It  Is  clear  that,  for  the  same 
mass  of  acid,  the  associated  solution  is  -harac cerized  by  a  greater  water  activity 
and  necessarily  water  vapor  pressure.  This  indicates  that,  when  the  primary  H^PO^ 
molecules  are  partially  associated,  the  vapor  pressure  lowering  effect  is 
diminished.  As  a  result,  the  water  vapor  pressure  gradients  between  solution 
droplet  and  environment  are  less  severe  for  the  condensed  phosphoric  acid  solutions, 
an  indication  of  reduced  hygroscopic ity . 

In  summary,  the  associated  thermodynamic  solution  properties  are 
established  from  the  unessociated  thermodynamics  through  molecule  association  and  by 
neglecting  the  chemical  differences  between  the  original  species  and  the  associated 

species. 


Water  Activity  of  Various  Condensed  Phosphoric  Acid  Solutions 
as  a  Function  of  H„PO,  Association 

j  'I 


Figure  3. 


3.2 


Equivalent  H,PO^  Phosphorus  Smoke  Acid  Concentration. 


By  definition,  the  a  id  i  mcentrat lor.  of  a  multicomponent  acid  mixture  Is 

given  as 


1+2,1, 31-Hni 
fMi+2,i,3i+ini  +  '"HjO 


(13) 


where  M^_2  ^  3i+l  1®  t^*e  molecular  weight  of  the  ith  order  condensed  phosphoric 
acid;  i  *  l’,  *being  H^PO^ ,  is  represented  by 

In  this  section,  the  equivalent  H3PO4  concentration  is  formulated  in 
order  to  compare  experiment  and  theory.  The  equivalent  ortho  acid  concentration 
function  assumes  that  each  mole  of  NaOH  dissolved  in  the  acid  solution  can  be 
equated  to  one  HjPO^  molecule  as  is  customary  in  titration  analysis. 

For  the  case  of  the  condensed  phosphoric  acids,  the  ith  order  condensed 
phosphoric  acid  will  possess  i  titration  end  points.7  Therefore,  the  total  number 
of  NaOH  molecules  dissolved  in  the  condensed  phosphoric  acid  solution  is 


[NaOH]  -  |inA  (14) 

Since  each  molecule  of  NaOH  is  equated  with  a  compound  having  a  molecular  weight  <.f 
H-jPO^,  the  equivalent  H-^PO^  mass  is 


Pi+2,i,3i+lni  *  M3 1 4? ln  i 


Substitution  of 
concentration  of 


equation  15  into  equation  13  results  in  the  equivalent 
an  arbitrary  mixture  of  condensed  phosphoric  acids  as 


1!,P0, 
5  “i 


F 

e 


a  M....U-X) 

~  3.14 


M  x  +  a  (1-x) 

H20*  -  314  * 


(If.) 


where  X  is  given  by  equation  4.  Figure  4  shows  the  dependence  of  equivalent  H-jI'o., 
concentration  on  relative  humidity  for  various  condensed  phosphoric  ;j;-M 
solutions.  As  the  solution  becomes  more  condensed,  F  increases  at  constant 
relative  humidity.  This  is  an  outgrowth  of  the  reduced  water  ahnorblnp 
characteristics  of  the  condensed  phosphoric  acids.  When  the  association  parameter  a 
is  equal  to  unity  (pure  H3PO4),  Ffi  reduces  to  equation  6  as  expected. 


1  7 


Figure  4.  Equivalent  H^PC^  Acid  Concentration  of  a 
Phosphorus  Smoke  Versus  Degree  of  Acid  Association 


Now,  for  the  present  case,  the  mole  ratio  of  pyrophosphorlc  acid  to  H3PO4 
is  expressed  in  terms  of  the  weight  percentages  as 


n 

n 


2 

1 


427 

31 


M_ 


314 


427 


07) 


where  and  are  the  weight  percentages  of  pyro  and  ortho,  respectively.  In 
putting  the  weight  percentages  experimentally  measured  by  Kaptshev,^  the  mole  ratio 
i  s 

—  -  1.65  (18) 
n314 

Now  each  pyrophosphorlc  acid  molecule  converts  to  two  molecules  through  branch 
point  hydration*^  The  corresponding  mole  ratio  of  acid  mixture  solute  number  2n^  to 
pure  hjP°4  solute  number  is 

ilni/?nt  ’  1,62  <>9) 
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But  this  is  exactly  the  association  parameter  a_  defining  the  reduction  in  solute 
number  due  to  the  association  of  H-jPO^  molecules.  Substitution  of  this  value  of  a_ 
into  equation  16  produces  the  equivalent  H^PO^  acid  concentration  of  a  phosphorus 
smoke.  Figure  1  shows  the  good  agreement  between  the  experimental  values  of 
Anderson  and  those  predicted  by  the  mixture  model  over  the  relative  humidity  range 
from  252  to  902. 

3 . 3  Real  Phosphorus  Smoke  Acid  Concentration. 

The  equivalent  H.PO^  acid  concentration  function  Fe  misrepresents  the 
real  acid  mixture  concentration  F_.  Although  the  equivalent  acid  concentration  does 
account  for  the  reduction  in  solute  hygroscopicity  through  solute  number  reduction 
in  the  case  of  theoretical  treatment  and  by  measuring  the  reduced  amount  of  water 
absorbed  in  the  case  of  experiment,  both  treatments  fail  in  describing  the  real  acid 
concentration  by  imposing  a  one  to  one  correspondence  between  NaOH  and  H^PO^.  If 
one  inspects  the  pyrophosphoric  acid  molecule  then  it  is  clear  that  two  NaOH 
molecules  are  associated  with  a  mass  which  is  less  than  two  molecules. 
Specifically,  the  pyrophosphoric  acid  molecule  is  represented  by  The 
result  is  that  the  equivalent  acid  concentation  function  overestimates  the  teal 
mass  of  acid  in  the  phosphorus  smoke.  More  mathematically 


fMi+2,i,3i+lni  M314llni  (20) 


The  real  acid  concentration  function  may  be  expressed  in  term6  of  the  equivalent 
ortho  acid  concentration  function  by  accounting  for  the  molecular  weight  differences 
between  multiples  of  H^PO^  molecules  and  pH  equivalent  condensed  phosphoric  acids. 

Beginning  with  the  equivalent  acid  concentration  Fe 


F 

e 


M3U  $lni/MT 


and  the  real  acid  concentration  function  Ff 


(21) 


Fr  "  5Mi+2,i,3i+lni/MT 


the  ratio  of  the  two  functions  can  be  expressed  as 


lMl+2,lt31+lni 

M31/,flnl 


(22) 


(23) 
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The  sum  in  the  numerator  of  equation  23  may  be  rewritten  from  equation  9  as 


IM  n 

i  i+2, i , 3i+l  i 


(M 


314 


-  M  ) 

h2o 


Ein  + 
i  1 


M  Er 
H„Oi 


(24) 


Combining  equations  23  and  12  results  in 


F 

r 

F 

e 


M 

113 


Sini  +  V*i 

M314  iini 


(25) 


Finally  substitution  of  equation  11  into  equation  25  gives  the  defining  relationship 

between  F  and  F  as 
r  e 


F  -  F 
r  e 


M  H  /a 

113+  H20  ~ 


M 


314 


(26) 


For  any  associated  solution  ( a^  >  1),  the  fraction  is  less  than  unity  and  the 
equivalent  acid  concentration  exceeds  the  real  acid  concentration.  Figure  5  depicts 
the  percent  deviation  between  Fg  and  Ff  for  a  series  of  condensed  phosphoric  acid 
solutions.  The  deviations  increase  monotonically  with  the  association  parameter^, 
varying  from  zero  at  a  *  1  to  -16%  at  £=  4.  Consequently,  the  greater  the  degree 
of  acid  association  the  greater  the  misrepresentation  of  the  real  acid 
concentration. 


(ASSOCIATION  PARAMETERS 

Figure  5.  Percent  Deviation  between  Equivalent  and  Real  Acid 

Concentration  as  a  Function  of  Acid  Association 


For  the  present  case,  _a  1.62  and 

F  -  0.93  F  (27) 

r  e 

Consequently,  procedures  involving  titrations  against  standard  0.01  N  NaOH  solutions 
will  overestimate  droplet  acid  concentrations  by  as  much  as  7%  when  a  one  to  one 
correspondence  between  moles  of  NaOH  and  H^PO^  ia  assumed. 


3.4 


Yield  Factor  of  a  PhosphoruB  Smoke 


Following  logic  similar  to  the  derivation  of  the  H^PO^  smoke 
(equation  7),  the  yield  factor  of  a  multicomponent  acid  smoke  is  given  by 

5Mi+2,i,3i+ini 

y.f.  -  - rr1— '-c -  (2b) 

*  Fr 

Here  the  phosphorus  aerosolization  parameter  z  and  the  composition  parameter  n  are 
set  to  unity.  Substituting  for  Fr  from  equation  23,  the  yield  factor  may  be 
expressed  as 


y.f. 


lin 
i  i 


314 


(29) 


Since  1  phosphorus  atoms  are  tied  up  in  the  ith  condensed  phosphoric  acid,  then 


My  -  M(P)  |in1 


(30) 


where  M(P)  is  the  mass  of  a  phosphorus  atom.  The  yield  factor  of  a  multicomponent 
phosphorus  smoke  becomes 


y.f.  =  M3U/M(P)  .  Fe 


(31) 


Comparison  of  equation  31  with  the  experimental  yield  factors  of  Tarnove^  indicates 
good  agreement  over  the  range  of  relative  humidities  from  40%  to  70%. 

In  conclusion,  it  haB  been  demonstrated  that  the  equilibrium  properties 
of  a  phosphorus  smoke  cannot  be  predicted  from  a  thermodynamic  model  which  is  based 
on  the  water  activity  of  H,P0,  solution.  Based  on  high  temperature  P,0,„  studies, 

J  H  4  1U 

chemical  composition  of  phosphorus  smoke  droplets  is  proposed  to  be  a  mixture  of 
orthophosphoric  and  pyrophosphoric  acid.  Applying  a  solute  number  reduction  scheme, 
the  solution  thermodynamics  of  the  condensed  phosphoric  acids  is  formulated  from  the 
ortho  water  activity  data.  Comparison  of  mixture  model  synthetic  data  and 
experimental  data  indicates  good  agreement. 
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4. 


CONCLUSIONS 


4.1  A  thermodynamic  model  which  Is  based  on  the  water  activity  of  H^PO^ 
solutions  will  underestimate  phosphorus  smoke  acid  concentrations  by  as  much  as  25% 
for  relative  humidities  less  than  90%.  Furthermore,  phosphorus  smoke  yield  factors 
are  overestimated  by  as  much  as  25%  for  relative  humidities  less  than  90%. 

4.2  Based  on  a  solute  number  reduction  scheme,  the  water  activity  of  ..n 
arbitrary  mixture  of  condensed  phosphoric  acids  is  formulated  from  existing  H^PO^ 
water  activity  data. 

4.3  High  temperature  P^O^q  hydration  studies  suggest  the  chemical  composition 
of  phosphorus  smoke  droplets  is  a  mixture  of  H^PO^  and  pyrophosphoric  acid, 
whereupon  the  solution  thermodynamics  of  the  acid  droplets  Is  developed. 

4.4  Comparison  of  mixture  model  synthetic  data  and  experimental  data 
indicates  the  condensed  acid  model  predicts  the  acid  concentrations  and  yield 
factors  of  a  phosphorus  smoke  within  5%. 

4.5  Due  to  the  molecular  weight  differences  between  the  condensed  phosphoric 
acids  and  the  pH  equivalent  multiple  of  H^PO^  molecules,  equivalent  H^PO^ 
concentrations  will  overestimate  the  acid  concentrations  of  a  phosphorus  smoke  by 

7%. 
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association  parameter 
water  activity 
droplet  diameter 

orthophosphor ic  acid  concentration 

equivalent  ortho  acid  concentration 

real  acid  concentration 

order  of  acid  association 

molecular  weight  of  orthophosphoric  acid 

molecular  weight  of  phosphorus 

mass  of  phosphorus  in  starting  material 

molecular  weight  of  water 

total  droplet  mass 

mass  of  water  in  droplet 

moles  of  water  in  droplet 

moles  of  the  ith  order  condensed  phosphoric  acid 
pressure 

droplet  saturation  ratio 
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ATTN:  DR  SMI  -  RPR  CDocuments)  1 

Reittone  Arsenal,  AL  35809 

L.  i  rector 

DAROUM  Field  Safety  Activity 

ATTN.  CK(4)S-u  I 

Ji  ar  I  oSTCwn,  IN  4  7  '  11 


Com  no  nder 

Ua  Army  Natlc*.  Research!  4 
Development  Qanriand 

ATTN:  DRDNA-VR  ' 

ATTN:  0RDNA-VT  I 

Nat  I cis ,  MA  0 1760 

US  ARMY  ARMAMENT  ,-ESEARCH  AND 
DEVELCFWENT  COMMAND 

Comrander 

US  Army  Armament  Research  and 
Dovelopnent  CttYitmnd 

ATTN;  0RDAR-LCA-1  1 

ATTN:  DREAR-ICE  ' 

ATTN:  DRDAK-LCE-C.  1 

ATTN:  OROAR-lCU  1 

ATTN:  DRDAR-LCU-CE  1 

ATTN:  DRDAR-FMA  (G.R.  Saccc)  i 

ATTN:  DRDAR-SCArW  1 

ATTN:  DROAR-TSS  2 

ATTN:  DftCPM-  CAWS- AM 
ATTN:  CR CPM-CAWS-S  I 

Cover,  NJ  07801 

Di  rector 

Ba  i  i  I  si  Ic  Research.  Laboratory 
ARR  ADuJd 

ATTN:  DRDAR-TSB-S 

Aberdeen  Proving  abound,  M)  21005 

US  ARMf  ARMAFi  NT  MATERi  EL  READ  I  NESS 

ccmmand 

Comnender 

US  Army  Armament  Materiel 
Fteadlness  Cbmnend 
ATTN-  DR  SAR- ASN 

ATTN:  CRSAR-PEM 
ATTN:  DR  S  ,AR-  S  F 

Hock  Island,  IL  61299 

Oomimnder 

US  Army  Dugway  Provl  ng  Ground 
ATTN:  Technical  Library 

Docu  Sect 
Dugway,  UT  8402  2 

US  ARMY  TRAINING  4  DOCTRINE  COMMAND 
Comnandant 

u5  Ai  my  iiifaritry  School 
ATTN:  NBC  Dlv I s Ion 

fort  tfennlng,  GA  31905 
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Commandant 

USAfP4CS/TC4FM 

ATTN:  AT2N-CM-CDM  1 

fort  McClel  Ian,  AL  36205 

Commander 

US  Army  Infantry  Center 

ATTN;  ATSH-CD-MS-C  ' 

fort  Efennlng,  GA  31905 

Commander 

US  Army  Infantry  Center 
Directorate  of  Plans  &  Training 
ATTN:  AT2B-DPT-PONBC  ’ 

fort  fanning,  GA  31905 

Commander 

USA  Training  and  Doctrine  Commend 

ATTN.  ATCO-Z 

Fort  Monroe,  VA  23651 

Commander 

USA  Combined  Arms  Center  and 
fort  Leaventorth 
ATTN:  ATZL-CA-OOG 

ATTN;  ATZL-CAM-IM 
Fort  Leavenworth,  KS  66027 

Commander 

US  Army  TRAOQC  Systen  Analysis 
Act  Ivi  ty 
ATTN:  ATAA-SL 

Vtiite  »nds  Missile  Ffcnge,  fW  68002 
US  ARW  TEST  4  EVALUATION  COMMAND 
Commander 

US  Army  Test  4  Evaluation  Command 
ATTN.  CRSTE-CM-F 
ATTN;  ORSTE-CT-T 

Aberdeen  Proving  Ground,  MD  21005 


Commander 

Naval  Weapons  Center 
ATTN:  Tecbn  Ice  I  Library 

Code  343  1 

CM  na  Lake,  CA  93555 

Commander  Of  f  leer 

Navel  WBapons  3jpor1  center 

ATTN:  Obde  5042 

Dr.  B.E.  Doudo  ' 

frane,  IN  47522 

US  MARINE  COFPS 

Director,  Development  Center 
Mar  I  ne  Corps  Deve  lepment  and 
Education  Command 

ATTN:  FI  re  Power  Dlvl  s  Ion  1 

Quantlco,  VA  22134 

DEPARTMENT  OF  THE  AIR  FORCE 

HO  foreign  Technology  Division  (AFSC) 

ATTN:  TQTR 

W-  I  jit -Patter son  aFB,  Oh  45433 
HQ  AFLC/LCWMM 

V4  i  jit -Patterson  AFB,  OH  45433 

OUTSIDE  AGEfCl  ES 

Battelle,  Co  lumbu  s  Labora  or  !  es 
ATTN:  TACTEC 
505  King  Avenue 
Qj  I  unbus,  OH  43201 

Toxicology  Inf  onrat  Ion  Centor, 

WAG  1 008 

National  Research  Council 
2101  destitution  Ave.  ,  Mm 
Washington,  DC  20418 


0ERAR1MENT  OF  THE  NAVY 
Comma  nder 

Naval  Explosive  Ordnance  Disposal 
foci  1 1 ty 

ATTN.  Army  Chemical  Officer 
d>de  AC-3 

Indian  Head,  kC  2C640 
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